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Abstract. The spectral theorem is proven for the quantum dynamics of quadratic 

open systems of n fermions described by the Lindblad equation. Invariant eigenspaces 
•n . of the many-body Liouvillean dynamics and their largest Jordan blocks are explicitly 

constructed for all eigenvalues. For eigenvalue zero we describe an algebraic procedure 
r^ ' for constructing (possibly higher dimensional) spaces of (degenerate) non-equilibrium 

steady states. 
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1. Introduction 

We study spectral decomposition of a general quantum Lindblad equation [H El [3] with 
quadratic generators [U |5] for an open system of n fermions. The physical picture and 
the motivation for the study of such a non-equilibrium quantum many-body problem 
have been outlined in Ref. |6] (see also j^), on which we shall rely extensively. However, 
in this note we show that an important assumption of diagonalizability made in Ref. 
[6] can be removed and that one can prove the spectral theorem for a general quadratic 
many-body Liouvillean and construct its canonical Jordan form. 

The main results presented bellow, are: (i) computationally efficienl||| canonical 
representation of the structure matrix of the Liouvillean quadratic form as summarized 
in lemmas 12. 4[ 12.51 (ii) writing the general (possibly non-diagonalizable) normal form 
of the quadratic Lioivilleans in terms of normal master mode maps [formula ( TMl) ]. (iii) 
constructing the full many-body Liouvillean spectrum, the invariant subspaces, and the 
largest Jordan block in each of them (theorem 14.11) . and (iv) algebraic construction of 
the convex subspaces of possibly degenerate (multiple) non-equlibrium steady states 
(theorem 14.20 . 

A general master equation governing time evolution of the density matrix p{t) of 
an open quantum system, preserving the trace and positivity of p, can be written in the 
Lindblad form [1] as 

^ = £p := -\[H, p] + Y^ {2L,pLl - {LlL„ p}) (1) 

where if is a Hermitian operator (Hamiltonian), [x,y] := xy — yx, {x,y} := xy + yx, 
and L^ are arbitrary operators, physically representing couplings to various baths. 

We shall treat open quantum systems of a finite number, say n fermions, with a 
quadratic generator £. Mathematically speaking, we shall restrict our discussion to the 
case where a set of 2n anti-commuting, Hermitian, bounded operators exist, say Wj, 
j = l,...,2n, 

{wj,Wk} = 25j^k, j,k = l,...,2n, (2) 

so that the Hamiltonian, and the Lindblad operators, can be simultaneously expressed 
in terms of a quadratic form, and linear forms, respectively 

2n 

H = y^ WjHjkWk = ui-Hw, (3) 

j,k=i 

2n 

X Computational problem of diagonalizing 4n x 4n com.plex matrix of Ref.[6i has been reduced to a 
problem of diagonalizing a,2nx2n real matrix (or possibly finding its Jordan canonical form if it cannot 
be diagonalized) and/or solving a 2n x 2n linear matrix equation of Sylvester type [S]. 
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Throughout this paper x = (xi, X2, ■ ■ .)^ will designate a vector (column) of appropriate 
scalar valued or operator valued symbols x^ and upright bold letter (e.g. H) will 
designate a complex or real matrix. 

The main conceptual tool of our analysis is a Fock space /C of operators (which 
may be referred to also as Liouville-Fock space) which the density operator p{t) ([I]) is 
a member of. /C is in fact a 2^" = 4" dimensional linear vector space spanned by a 
canonical basis |Pa) with 

Pa„a„...,a,^ ■= T'^'^wTw'i'^ ■ ■ ■ w^" , a, E {0, 1}, (5) 

which is orthonormal with respect to Hilbert-Schmidt inner product 

{x\y) = iix'^y. (6) 

Introducing creation and annihilation linear maps Cj over /C 

c]\Pg_) = 5a^fi\WjPa) Cj\Pa) = 5a^^l\WjPa) (7) 

the space /C becomes a fermionic Fock space. Canonical anticommutation relations 
(CAR) among operator maps {cj, c^} = 0, {cj, c|.} = 5j^k and the fact that c] is indeed 
a Hermitian adjoint of Cj w.r.t. (|6]) can be checked straightforwardly. 

As shown in Ref.jB] the Lindblad equation ([T]) in the Liouville-Fock picture 

^|P) = ^|P) (8) 



takes manifestly quadratic form, where the Liouvillean map is expressed as | 

£ = a- Aa-Aoi (9) 

in terms of An Hermitian Majorana fermionic maps over the operator space /C, 

aij := {cj + c])/V2, aa,, := i(c,- - c])/y2 (10) 

satisfying CAR 

{O'u^j.af.^k] = 5u,t,5j,k^ z/,yu = l,2, j,k = l,...,2n. (11) 

We note that physically interesting non- equilibrium steady state (NESS), |NESS) e /C, 
can be defined as the (possibly degenerate, and then non-unique) right zero-eigenvalue 
eigenvector of the Liouvillean £|NESS) = 0. The antisymmetric 4n x 4n structure 
matrix A ([9]) has the general form [6j 

-2iH + 2iMi 2iM 

-2iM^ -2iH-2iMi 

= - 2il2 ® H - 2(7^ ® Mr - 2((T^ - ia^) ® M; (12) 

where a", a = 1, 2, 3 are the Pauli matrices, H = — H^ is an antisymmetric (due to (|2])) 
2n X 2n Hermitean (purely imaginary) matrix, whereas 

Mr := i(M + M) = M^, Mi := 4(M-M) = -Mf (13) 

§ Note also a technical issue of a parity symmetry of the Liouvillean and the operator space w.r.t. even 
and odd number of fermionic operator-excitations W which we may ignore here without affecting any 
of the results. 
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are the real and imaginary parts of the positive semidefinite Hermitian 2n x 2n Lindblad 
hath matrix M := 'Y^ in ® I-, M > 0. Note that, physically speaking, H encodes the 
Hamiltonian, Mr encodes the dissipation and Mi encodes the driving. The scalar Aq 
(|9]) has a value 

Ao = 2trM = 2trMr. (14) 

Writing the skew unit 

J := a^ ® l2n (15) 

we observe that the structure matrix satisfies a particular self-conjugation property 

A = JAJ, or equivalently AJ = JA. (16) 

2. Preliminaries 

A series of useful lemmas on the spectral decomposition of the structure matrix A shall 
now be straightforwardly proven, which explore the structure of the problem and go 
beyond the general spectral theorems for the anti-symmetric complex matrices [9]. 

Lemma 2.1 Structure matrix [W\) is unitarily equivalent to a block-triangular matrix 

A = -2il2 (g) H - 2a^ ® Mr + 4i(T+ O M; = f ~^ ^^^' j (17) 

where 

X:=-2iH + 2Mr (18) 

is a real 2n x 2n matrix. 

Proof The similarity transformation A — t- A 

A = UAU^ (19) 

is trivially provided by the cyclic permutation a^ — ;■ cr^, cr^ — )■ cr^, a'^ — )■ cr^. 

We shall refer to the representation in which the structure matrix has a form (TT7|1 
as the tilde-representation. As an immediate consequence of lemma 12.11 we have: 

Corollary 2.2 The characteristic polynomial of the structure matrix A factorizes in 
terms of the characteristic polynomial of X fT8\} 

PA{f3)=p^{f3)p^hf3), (21) 

where Pa{P) '■= det(A — /3l4„),px(/3) := det(X — /3l2n) o^nd hence, if a (A) denotes the 
spectrum of a matrix, 

a(A) = a(X) U (-a(X)), (22) 

where the union U should keep track of (algebraic) multiplicities. 
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Curiously enough, this imphes that the spectrum of the structure matrix is insensitive to 
the driving matrix Mi. Note also that the real matrix X flTSj) has no general structure 
apart from the fact that X + X^ = 4M,. > 0, since M^ > is implied by M > 0. 
Sometimes the eigenvalues {/3j} of X shall be referred to as rapidities. 

Lemma 2.3 Let ^ be a real square matrix, such that X + X-^ > 0. Then: 

(i) Any eigenvalue (3 ofX. satisfies Re/3 > 0. 

(a) For any eigenvalue (3 of X. on the imaginary line, Re/3 = 0, its algebraic and 
geometric multiplicities coincide. 

Proof (i) Take an eigenvalue /3 and an eigenvector u, write Xm = /3m, and the complex 
conjugate of this equation Xm = /3m. Then take a dot product of the first equation with 
u and the dot product of the second equation with u and sum up: 

m-(X + X^)m= (2Re/3)M-M. (23) 

Strict positivity of the eigenvector norm, u-u>0, and non-negativity, u- (X + X"^)m > 
(by definition of positive semidefiniteness of X + X^), imply Re/3 > 0. 
(ii) Consider a linear system of differential equations, 

{d/dt)u{t) = -Xu{t). (24) 

Positive semidefiniteness of X + X^ is then equivalent to Lyapunov stability in control 
theory, namely (d/dt)||M||| = —u ■ (X + X-^)m < iff X + X"^ > 0. Then, the converse 
(equivalent) version of the statement (ii) of the above lemma says that: If there exists 
an imaginary (or vanishing) eigenvalue /3, /3 = ib, and a corresponding Jordan block 
of dimension A; > 1 in the Jordan canonical form of X, then X + X"^ is not positive 
semidefinite. This is indeed obvious, since if we take the initial vector m(0) for fl24|) from 
ker (X - (3lf Q ker (X - f3lf-^ then u{t) oc t'^-^e"'^*, and so X + X^ ^ 0. | 

Lemma 2.4 Let us write the Jordan canonical form of the matrix X ( fl^) as 

X = PAP-^ (25) 

where P is a non-singular 2n x 2n matrix (columns of which consist of generalized 
eigenvectors ofX.) and A = ^-^ A^^ ,,(/3j) is a direct sum of Jordan blocks 

//3 1 \ 



MP) ■-- 



/3 ■■■ 
■■■ 1 

V /3y 



(26) 



for the distinct eigenvalues - rapidities f3j, and the block sizes ij^k satisfying the 
completeness sum 






; ^ ij,k = 2n. (27) 
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Then, the manifestly antisymmetric Jordan canonical form of the structure matrix [W^] . 
specialized to the case of zero driving Aq := A|Mi=oj 

/ -2iH 2iMr 

Ao= ( _ 

reads 



-2iMr -2iH ' ' *^^^ 



Ao = Vo^ ( _^T ) Vo (29) 



where 



Vo = (P^ © P-')U (30) 

is a An X An matrix (rows of which store the generalized eigenvectors of Aq) which 
satisfies normalization identity 

VoV^ = J. (31) 

Proof Specializing lemma [2711 to the case of Mj = 0, we have Aq = (— X"^) ©X. Using 
the Jordan decomposition fl25|) and the similarity transformation flT9|) we find 

Ao = Ut(P-^ © P)[(-A)^ © A](P^ © P-^)U. (32) 

Introducing Vq := (P"'" © P^^)U, i.e. equation (!30|) . and noting Ut = U"^J, we write 

V^i = Ut(P~^ © P) = U^(P © P-^)J = V^J (33) 

which immediately implies ( l29l) and (pTl) . | 

Lemma 2.5 The general structure matrix A 17^) is always similar to its zero-driving 
counterpart Aq ( f^) 

A = W-^AqW. (34) 

The similarity transformation is provided by a complex orthogonal matrix 

W = l4„, + 2((T^ - ia^) © Z, WW^ = l4„, (35) 

where Z is an antisymmetric 2n x 2n real matrix, TF = —Z, obtained as a solution of 
the continuous Lyapunov equation [lO^ 

X^Z + ZX = Mi (36) 

which always exists (for the straightforward method of solution, see the proof). 

Proof Working in the tilde-representation ( 1T^ . W := UWU''' = l^n — Aia'^ © Z, we 
have W"^ = 1^^ + 4icr+ © Z and 
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so fl34p is indeed consistent with f ll2p iff Z is antisymmetric and solves the continous 
Lyapunov equation fl36|) . [jj] 

We shall now show by construction that under condition of positive semidefiniteness 
of M = Mj. + iMi, the solution of (136!) always exists, but may not be unique if some of 
the eigenvalues /3j of X (IT8l) lie on the imaginary line. 

Writing a trivial mapping vec : M2nx2n _^ ]^4n ^ which vectorizes a matrix 
(vecX)2„(j_i)+j := Xjj, and obeys vec (XZY) = (X (g) Y-^)vec Z, the continuous 
Lyapunov equation ( l36l) reads 

(X^ (g) l2„ + l2„ ® X^) vec Z = vec Mj. (39) 

This equation is solved efficiently (in 0{rc') steps) using the Jordan canonical form of 
X ([25]) 

( A^ ® l2„, + l2„ ® A^) vec (P^ZP) = vec (P^MjP) (40) 

Note that the (2n)^ x (2n)^ matrix f2 := A-'" ® l2n + l2ra ® A-^ is lower triangular with 
diagonal entries fi,. ,, of the form /3j + /3j/, so vec (P-^ZP) can be solved for uniquely if no 
such pairs /5j, /3j/ exist that /3j + Pji = [HI [11]. Since Re/3j > the problem {flr,r = 0) 
may arise in either of the two cases: (i) (3j = 0,(3j' = 0, or (ii) /3j = ib,(3j' = —ib, 
b e M\{0}. It follows from lemma [231 that in either of these cases (i,ii) the corresponding 
Jordan blocks are trivial (of dimension 1), so there are no non- vanishing elements in the 
entire r-th row or column, Q^.s = ^s,r = 0. Thus there will be a solution of ([40]) only if 
the corresponding coefficient on the RHS vanishes 

a;, :={vec(P^MiP)}^ = 0, (41) 

and then, the solution may be non-unique as the coefficient {vec (P-^ZP)},, is arbitrary 
(unless it is fixed to zero by the antisymmetry of Z like in the case j = j', (3j = 0). Let 
us now study these two singular cases: 

(i) (3j = and /3j> = 0. Suppose eigenvalue has multiplicity m = dimkerX, 
and denote the corresponding basis vectors Ae^ = 0, ase^, I = l,...,m. The 
RHS coefficient Ur pi]) can be any of the matrix elements of an antisymmetric 
m X m matrix Ki^i, := (e^ (g) e^,) ■ vec (P'^MjP) = (Pe^) ■ Mi(Pe;,) = -i^r,«- 
Note that the corresponding eigenvectors of X, Pe^, are real (since X is a real 
matrix), and hence the matrix K is real as well. Positive semidefiniteness of 
M implies positive semidefiniteness of its counterpart reduced on the null space 
Qi^v := (Pe,) ■ M(Pe,) = (Pe,) ■ {i(X + X^) + iMi}(Pe,) = ii^^', Q > 0. Since 
K"^ = — K, this is only possible if Kii' = 0. 

II All admisible structure matrices (J12p form a non-semisimple n(6n — 1) dimensional Lie algebra 
spanned - in a tilde-representation - by 4n x An matrices of the form 

-X^ iY' 



X J (^«) 

where X is an arbitrary 2n x 2n real matrix and Y is an antisymmetric real matrix Y — — Y. The 
similarity transformation (IMI) is generated by a nilpotent subalgebra with X = 0. 
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(ii) Pj = ib, f3j/ = —16, 6 G M \ {0}. Suppose eigenvalue ib has multiplicity m = 
dimker (X ± ibl2n), and denote the corresponding basis vectors A/ = ±16/ , as 
/ , / = l,...,m. Since X is real, the corresponding eigenvectors uf" '■= P/, , 
^uf" = ±i6wf, can be chosen to be conjugate paired u^ = y^ (i.e. the 
corresponding columns of P can be chosen to be conjugate paired). The RHS 
coefficient cu^ (jH]) can be any of the matrix elements of an antihermitian m x m 
matrix Ki^i> := (/+ ® /p ■ vec (P^MiP) = u/ ■ Mim^T = -Ki',i. Positive 
semidefiniteness of M^ := M and M~ := M implies positive semidefiniteness of two 
Hermitian mxm matrices Q^, := m/^-M^m^T = M/'-{|(X+X-^)±iMi}M^7 = ±17^^^^/, 
Q''^ > 0. iK and — iK can be simultaneously positive semideifinite only if Kii/ = 0. 

I 

3. Jordan canonical form on the tensor product spaces 

To best of the author's knowledge the theory of Jordan canonical form on tensor 
product spaces has not been very much developed. However, being able to manipulate 
Jordan decompositions on tensor product spaces seems crucial for discussing the general 
non-diagonalizable situation of quantum Liouvillean dynamics. In this section we are 
demonstrating two simple results, namely one on Jordan decomposition of a sum of 
Jordan blocks on tensor product spaces, and the other one on Jordan decomposition of 
a many-body nilpotent map on the fermionic Fock space whose structure matrix is a 
simple Jordan block. Both results will be needed for the proof of the spectral theorem 
later in section HI 

3.1. Jordan canonical form of the sum of Jordan blocks on a tensor product space 

Let us consider a tensor product space C*^ C?) C' and ask for the Jordan canonical form 
of the following kl x kl matrix 

Ak,i = Ak{a) ^ h + Ik ^ Ai{f3) (42) 

where the Jordan block matrix is defined in (l26l). 



Lemma 3.1 The matrix [4^ has the following Jordan canonical form 

rmn{k,l} 

Ak,i- Afc+,_2r+i(« + /3). (43) 

r=l 

The proof will be constructive, so we shall also show how to compute the corresponding 
generalized eigenvectors. 

Proof Without loss of generality we can assume k > I. Let \i,i),i = l,...,k,j = 
1,...,/ denote the basis of C'^ C^ Clearly, since Ar(7) = Ir + tA^ (we write 
Ar := Aj,(0)), it is enough to show (H5]) for the nilpotent case a = f3 = 0, and then at 
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the end add {a + /3)lki. The crucial tool will be the Newton's binomial formula on the 
tensor product sum 

r=0 

From the above formula, and nilpotency A![; = 0, Aj;~^ ^ 0, it follows that A^| ^^ = 
but A^y^^ = C'j^'^^) A^"^(8)A|"^ 7^ 0, so the /ar^esUordan block in the decomposition 
of Afc^i has size k + I — 1, with the first nilpotent series of generalized eigenvectors 

b;l):=A^JO;l), p = 0, . . . ,k + I - 2, (45) 

and a 'seed' vector (null vector of A^|'^^), |0; 1) := \k,l). 

Let us introduce the diagonal number matrix N := J2i=i J2j=ii'^~^J)\h J) {h j\ which 
slices the vector space to a direct sum of eigenspaces of the following dimensions 

dimker(N - (g + 1)1m) = dimker(N - (fc + / + 1 - q)lki) = q, (46) 

which change exactly by 1 when changing the eigenvalue by 1. Further, note that 
NAfc/ = Ak^i(N — Iki), or equivalently A^; always decreases the eigenvalue of N 
exactly by 1 

Afc,,ker(N - (g + 1)1^,) C ker(N - ql^i). (47) 

Then we make the following inductive argument. Suppose that for all q < r, the 
series of generalized eigenvectors 

b;g):=A^JO;g), p = 0, . . . ,k + I - 2q, (48) 

explore all the eigenspaces of N with eigenvalues < r and > k + l + 2 — r. This is trivially 
true for r = 2, or the maximal series with g = 1 which we have constructed above ( H5l) . 
We need to show now that the r-th series explores the other two extremal 
eigenvalues, r + 1 and k + I + 1 — r and is of length k + I — 2r + 1. Clearly, the 
seed vector for the r— th series |0; r) should have the number eigenvalue k + I + 1 — r, 
i.e. it should be of the form 

r 

|0;r) = ^c(")|fc-g + l,/-r + g), (49) 

9=1 

and the size of the corresponding Jordan block can not be longer than k + I — 2r + 1 

Ky-''->, r) = g |r - ,,,) ±(Y [- '; ,\') <^ = 0. (50) 

This results, for the rth Jordan block, in a homogeneous system of r — 1 equations for 



r unknown coefficients c^ 



(r) 



E (*'!'+ J) 4" = 0. ie{i....,r-i}. (51) 
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where on RHS and below we write k' := k — r, V := I — r. The above system has (up to 
a prefactor) a unique solution 



4"=(-ir' ::':': . (52) 



4" -Mr- (*'::;') ('';!;' 

What is left to show is that a nilpotent series starting from a seed |0;r) is indeed not 
shorter than fc + / — 2r + 1, i.e. A^"J"'~^''|0; r) ^ 0, or equivalently 

i:(;:;:>r^o. ,53, 

for at least some j G {l,...,r}. Verification of (15T!) and (153!) with ansatz (l52l) is an easy 
case for a computer-assisted proving system based on Zeilberger algorithm [12]. | 

5.^. Jordan canonical form of the many-body nilpotent map 

Here we study the following abstract many body problem. Let bk, b'j^,k = 1, . . . ,£ be 
fermionic maps satisfying CAR, {bj,bk} = {b'j,b'i.} = 0,{bj,b'f,} = 6j^k, and defining i 
fermionic modes, i.e. spanning a 2^ dimensional vector space V with the Fock basis and 
its dual defined, respectively, as 

\u,,...,u,):=ib\)'''---ib',f\nght), (54) 

(i/i, . . . , z/,| := (leftKfe,)'^^ ■ ■ ■ ihr\ ^k e {0, 1}. (55) 

The vector |right) and its dual (left|, (left|right) = 1, are defined via (left|6';. = 
0,6fc|right) = 0, hence the biorthogonality (i^|i^') = 5^,]/ simply follows from the CAR. 
We shall study the following nilpotent linear operator over V 

e.-i 

M = J2b'k+ibk. (56) 

fe=i 

We observe immediately that A4 commutes with a number map J\ := 'Ylik=ib'kbki 

[M.,J\f] = 0, which foliates the linear space to a direct sum V = 0^=0 "'^"i' ^^ 

dim Vm = (,^) = rn\(e-mV. dimeusioual spaces V„, which are spanned by Fock vectors {u) 

with exactly \u\ := J2k=i ^k = fn 'particle' excitations (m being the eigenvalue of A^). 

We shall now fix the number of particles m and restrict the nilpotent operator 
Ai f l56|) to the space Vm, M.m '■= A^|v,„- For m = (and m = £) A^ acts trivially, 
M.Q = Mm = 0, and for m = 1 (and m = £ — 1), Mi ~ Mi-i ^ A^. However, the 
interesting question now is to find a Jordan canonical form of Mm, for arbitrary m. 

Let us define a weight of a basis Fock vector \v) as 

^ 1 

w^= ^A;z/fc - -m(m + 1). (57) 

k=l 

Having (!.,m fixed, there are two extremal states, with minimum and maximum weight, 
namely 

|min) = |1, . . . , 1, 0, . . . , 0), Wmin = 0, 

|max) = |0,. . .,0, 1,. . .,1), Ujna^ = {i - m)m. (58) 
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Denoting the linear spans V^ = L{\iy); \u\ = m,u^ = r}, we have a direct sum 

{i—m)Tn 

V„.= V:^. (59) 

r=0 

Dimensionahty dimV^ =: (^) has a clear combinatorial-physical meaning, namely 
(^) is the number of distributions (combinations) of m particles on i exclusive 
equidistant vertical levels, such that the total potential energy is exactly r steps 
larger than the minimal. Such a restricted binomial symbol can be calculated from 
a straightforward recursion relation 

-1\ f i-l\ /O 



II I ■ \ 1 / > \ I ^m,oSr,0 (60) 

/ r \ / r \ / r—l-\-m \ / 

and clearly satisfies the sum rule 

' * \ m I \ m I 



r=0 \ / r \ 



and the symmetry 



' r \ / [l—m)m—r 



(62) 



Note that (,^) =0ifr<0orr> m{i — m). Crucial observation for the rest of our 
discussion is the fact that Aim maps V^ to V^^, i.e. that the map A4 fl56|) increases 
the weight of the state exactly by 1, 

MV:^ C v;+\ (63) 

Clearly, one can now state the following: 

Lemma 3.2 The largest Jordan block in the canonical form of Aim has size {i—m)m+l. 

Proof Taking |0; 1) := |min) as a 'seed' for the first nilpotent series, we see that, due 

to dnsD, 

|p; 1) := MP\mm) E V^. (64) 

All these vectors for p < {m — i)m should be non-vanishing since only positive integer 
coefficients can appear in expressing the series A^P|min) in the canonical basis lu), and 
the only canonical basis vector |z/) of Vm. which is annihilated by A4 is |max). So we 
have \{i — m)m; 1) = g|max), where g is a positive integer since dim Vm = 1, and 

M\{i-m)m]l) = 0. I 

The rest of Jordan decomposition of A^„ we are only able to state if the following 
conjecture is true: 

Conjecture 3.1 The map M. : V^ — )■ V^"^ is injective forr < [_{{£ — m)m — l)/2j and 
surjective for r > \_{i — m)m/2\ . 
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Conjecture has been verified explicitly using computer algebra code for all i < 12, and for 
any £ and m < 3 [13] , however its general proof remains open due to a rather involved 
combinatorics. The conjecture of course also implies the ordering on the restricted 
binomial sequence (dimensions of V^), namely 

n<(n<...<(n (65, 

which can be verified directly. 

Omitting a straightforward proof [I3] we state that the conjecture implies the 
following simple general result 

[(e-rn)rn/2\ 
-^" - 95 ~ ^{e-m)m+l-2r (66) 

where the notation [k]A. means A © A ... © A(A; times), and [0] A := 0. In other words, 
the number of Jordan blocks of size {i — m)m + 1 — 2r (it can only decrease from the 
maximal in steps of 2) is (^) — (^) _ (which can also be zero for some instances r). 
The total number of Jordan blocks and thus the total number of proper eigenvectors of 
Mm is then simply equal to (^) l(,_„)„/2j • 

4. The spectral theorem 

With the lemmas 12.4112.51 we have shown that any Liouvillean structure matrix of the 
form f lT2|) admits a canonical representation in terms of the Jordan canonical form 
( l25|26l) of the real matrix X ( IT8l) . and the anti-symmetric real matrix Z = —Z'^ which 
solves the Lyapunov equation (!36|) . namely 

A 

A^ 

The eigenvector matrix (rows of which store the generalized eigenvectors of A) reads 

V = VqW = (P^ © P-^)U(14„ + 2{(x' - ia^) © Z) 
1 /P^(l2n-4iZ) -iP^(l2„-4iZ)\ 

and satisfies the canonical normalization condition 

VV^ = J. (69) 

Let us name the first 2n rows of V as Vj^k,i^ ^^^ the last 2n rows as v^j^k,h which are 
exactly the generalized eigenvectors pertaining to k-th Jordan block of the eigenvalue 
(rapidity) f3j, and —(3j respectively, and / = !,...,£ j^k (^ = 1 designates the proper 
eigenvector) where ij^k is the size of the Jordan block {j,k). Then we introduce the 
normal master mode (NMM) maps as 

bj,k,i ■= Vj^k,i ■ a, b'j^k,i ■= v!j,k,i ■ «, (70) 

satisfying the almost-CAR (simply following from ( 169|) ) 

{bj,k,i,bj'^k',i'} = 0, {bj,k,i,b'j,k',i'} = Sj,j'h,k'Si,i', {b'j^kpb'j,k',v] = 0. (71) 



A = V^( ^^ , )V. (67) 
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In terms of decomposition f l67p and NMM maps (170]1 the Liouvillean Q gets almost- 
diagonal form 

C = (Va) • [A © (-A^)] J(Va) - A^i (72) 

j,k 1=1 

+ Yl ihk,ib'j,k,i+i - ^i,fc,z+i^i,fc,0} (73) 

i=\ 

where the scalar (TT^ has been expressed as Aq = trX = Xljfe^i.fc/^j- After exercising 
f lTTj) the Liouvillean finally simplifies to a normal form 

(3^J2^b'^Aikk,i+ E WA>4- (74) 

z=i z=i J 

Note that the second (non-diagonal) sum vanishes for the trivial Jordan blocks with 

kk = 1- 

Following Ref. [6] a unique paiiiJll of operator Fock states exist, satisfying 

6,-,,,|NESS) = 0, (1|6;.,,, = 0. (75) 

The bra (dual) vector (1| in fact represents a trivial operator 1, dual to |NESS), 
(1|NESS) = trpNESs = 1- Indeed |NESS), which is an empty Fock state w.r.t. 
creation NMM maps b'j^k,h i^ a non- equilibrium steady state (NESS) since straightforward 
inspection following from almost-diagonal form of C fl71|) shows that £|NESS) = 0, as 
well as the dual relation (1|£ = 0. 

The complete spectrum of eigenvalues and eigenvectors of 4" x 4" matrix of C over 
/C is trivial to construct in case that all Jordan blocks of X are trivial [ij^k = 1), i-e. 
if A is diagonalizable, the case which has been assumed in Ref. [6]. However, here we 
treat the general case. 

4-1- Spectral decomposition and Jordan canonical form of the Liouvillean 

Theorem 4.1 (i) The complete spectrum of Liouvillean C is given by the following 
integer linear combinations 

A„ = -2 E rrij^kf^j , rrij^k G {0, 1, . . . , ij^k}- (76) 

j,k 

(a) The 4" dimensional operator space, and its dual (the bra-space), admit the following 
decomposition 

/C = 0/C^, /C' = 0C (77) 



m 



% Unique with respect to a chosen, possibly non- unique (see lemma [53]) decomposition of the structure 
matrix (l67l). 
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in terms o/ dim /Cm = Y\.jk ( m'' ) dimensional invariant subspaces 

ClCrn ^ A^m , ^'rrA ^ ^m (78) 

spanned by 

n n ^;-,fc,,jNESS); l<h<...< ^, < £,.fc 

j,fc ri=l 

(ii n n ^^•-'^.'.5 1 < /i < . . . < /m,,, < £,,. . (79) 

j,k ri=l J 

(Hi) However, the dimension of the eigenspace (the number of proper eigenvectors 
corresponding to Xm) is smaller than dim/Cm in the nontrivial case when at least one 
ij^k > 1- The size of the largest Jordan block corresponding to Xm is 

1 + ^{ij,k - mj^k)mj^k- (80) 

Proof We start by showing (ii). The fact that direct sums ( 177|) span the entire space 
follows from dimension counting ^^ dim /Cm = X^mriiA; ( ''■'" ) ~ Ilifc^^^'* = 2^" and 
non-singularity of the complex orthogonal transformation to NMM maps (!70l) . The 
invariance conditions (1781) are equivalent to 

(y\L\x) = 0, V|x) G /Cm, {y\ € /Cm'; such that mj^rnf, (81) 

which in turn follows from the normal form ( 17^ . More explicitly, one may write fl74l) 
as £ = J2j,k ^j,k where 

4fc := -2{P,Arj^k + M,,k), (82) 

^j,k 3-,k -'- 

■^j^k = 2^ b'j,k,fijAh Mj^k = 2^ b'j,k,i+ibj,k,i 
1=1 1=1 

and show (IHTl) separately for each term {y\Cj^k\x) = 0, and considering only \x) = 
U7=i ^^-.fe.^jNESS), {y\ = (1| ll7=t kk,i, for some {/J, {/;}, since the other NMM maps 
bj',k',i' for {j',k') 7^ {j,k) entering {x\ and \y) (I79p all annihilate as they commute with 
Cj^k- The rest is just trivial CAR algebra. 

To show (i) we write the Liouvillean as £ = £o + -^ where £o := ~'^J2jkf^J-^i<k 
and Ai = J2jk-^j,k- Since Aij^k are nilpotent (lemma \^72\i . and [Aij^k,-Mj'^k'] = 0, 7W 
is nilpotent as well. Clearly [Cq, J\4] = as all the terms commute [Afj^k, ■^j',k'] = 0, so £ 
and £o should have identical spectra. But the eigenspaces of Cq are /Cm, £o/Cm = XmK-rn, 
so {Am} is also the spectrum of £. 

As for (iii) we observe that the invariant space /Cm admits a tensor product 
decomposition 

/Cm ~ Q^ /Cj,fc,mj^fc (83) 
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where lCj^k,m, is a {^l^) dimensional space spanned by all Y\^=i^'jki |NESS), with 
1 < Zi < . . . < /m < £j fc. Let us consider the restriction of the map on the invariant 
space jC|/c,„ = A^l + M.^ where the nilpotent M.m acts independently on tensor factors 




M.^=y:\ (g) U,,,^^,\®M^;i. (84) 

i'.fe' ' ' / 

A^m/fc is exactly the abstract nilpotent many-body map studied in subsection 13.21 where 
^ = ^j,k- Noting the lemma IX^ each factor map has a nontrivial Jordan canonical form 
with a maximal block of size {i.j^k ~ fT^j,k)^T^j,k + 1- Then we use lemma |3TT] recursively, 
noting at each step always only the largest Jordan block, to arrive at ( IHOl) . | 

Writing down the complete Jordan canonical form of Liouvillean C seems at this 
point a very difficult combinatorial problem even assuming the conjecture 13.11 were 
true. We stress however an interesting non-trivial aspect of the result that we 
have just proven. Namely, in spite of the fact that the dimension of the invariant 
subspace pertaining to some possibly highly degenerate rapidity /3j is exponentially 
large in the (algebraic) multiplicity of the rapidity, the largest Jordan block is only 
polynomial (at most quadratic) in the multiplicity. This may have interesting physical 
consequences in studying Liouvillean time evolution (relaxation) where all terms of the 
form t^"^ exp(Am_^) can appear with integer power p not larger than ( IHOl) . 

4-2. Uniqueness of NESS 

Theorem 4.2 Uniqueness of NESS. |NESS) is a unique stationary state of open 
quantum dynamics ^ if and only if all eigenvalues f3j of X ( fl^) lie away from the 
imaginary line Re/3j > 0. If this is not the case, then: 

(i) For each zero rapidity (3j = 0, 

|NESS;j,A;) :=6;.;._i|NESS) (85) 

we also have the stationarity £|NESS; j, k) = 0. 
(a) For each imaginary rapidity Pj =ib, b eW\ {0}, we have a corresponding negative 
rapidity (3ji = —ib, and 

|NESS;j,/, A;,fc') := 6;.,,i6;.,,,, JNESS), (86) 

which satisfies stationarity £ | NESS ;j,j', k, k') = 0. 

Proof Uniqueness of |NESS) in case of strict positivity of the spectral gap of £, 
A = 2minj Re/3j (theorem 14. ip is obvious. In the other cases, (i), the stationarity 
£|NESS; j, k) = follows from the facts that the term j, k is absent in the normal form of 
the Liouvillean (17il) . since /3j ^ = and £j^k = 1 (lemma [273]) . whereas b'^,^^^ anticommutes 
with all the other terms,. So the annihilation maps b can be commuted to the right 
and then ( 175|) is used. Similarly, in case (ii), applying C ( 17il) to \NESS;j,j',k,k') and 
commuting annihilation maps b to the right results in two more terms w.r.t. case (i), 
which however cancel each other since (3j + /3,/ = 0. I 
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Note that (1^ implies 

b = p-^c^ c^ = ^(fli + ia^) (87) 

where a^^, 03 represent, respectively, the first and the last 2n components of a 4n vector a. 
Therefore, since cj can be interpreted as creation maps of Hermitian fermionic Majorana 
operators Wj [6j, we have the following obvious statements: in case (i), the row (j, /c, 1) of 
P^^ which is the left eigenvector of a real matrix X of eigenvalue can be always chosen 
real, and hence |NESS; j, /c) represents a Hermitian operator , and similarly in case (ii), 
the rows (j, k, 1) and (j', k, 1) of P^^ can be chosen to be mutually complex conjugate 
(as they correspond to left eigenvectors w.r.t. a pair of mutually conjugated eigenvalues) 
and hence |NESS; j, j', fc, /c) represents a Hermitian operator. In other words, h'-^,^ (i) 
and b'jj^^^b'j,).^^ (ii) are, under the above conditions on the eigenvectors of P, Hermiticity 
preserving maps. Furthermore, using similar arguments one can see that the non- 
diagonal maps b'j^k,ib'j',k',i + ^j,fcM^i',fc,i and i{b'j^k,ib'j',k',i - ^i,feM^i',fe,i) ^^^ Hermiticity 
preserving, so 

\NESS; J, j',k,k',+):=mSS; J, f,k,k') + |NESS; j,/. A;', A;) 

INESS; J, /, k, k\ -) := i|NESS; j, j', k, k') - i|NESS; j, j'. A;', k) (88) 

correspond to Hermitian operators. 

Note as well that all these elements of /C correspond to trace-zero operators as 
we have trivially trpNESS;j,fc = (1|NESS; j, fc) = (1|6^-;, ^jNESS) = 0, trpNESS;jj',fc,fc' = 
(1|NESS; j, j',fc,fc') = (l|6^.;i^i6^.,^;.,_i INESS) = 0, due to (I75]). Thus, in general a convex 
set of real coefficients should exist {aj,fc, C(j,j',k,k',± £ ^}, such that 

INESS) + J2 «i,fc|NESS;j, A;) + ^ a,-,,,fc,,v|NESS; j,/. A;, A;', r) (89) 

j,k j,j',k,k',T=± 

represent valid positive trace-one density operators. 

Geometric characterization of this set of degenerate non-equilibrium steady states 
poses an interesting open problem. 

4-3. Physical observables in NESS 

Theorem 4.3 Suppose the spectrum of~K. ( I7g|) is strictly away from the imaginary axis 
Re/3j > 0. Then the complete set of quadratic correlators in NESS is given by the 
matrix Z, a (unique) solution of the Lyapunov equation ^3R) . 

trWjtiJfcpNESs = ^j,k + 4iZj,fc. (90) 

Proof Inverting ( 1701) 

a = Y~\b®y)=Y^ib®b), (91) 

we express the first first 2n components of the canonical Majorana maps explicitely 

ai = -^ JP-^ + (I2. + 4iZ)P6'} . (92) 
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Then we note that, since \wjWk) = 2dijai^k\^) and {wjWk\ = 2{l\aijai^k, we can express 
the correlation matrix Cj^k := trWjWkPNESS as 

C = 2(l|ai®ai|NESS) 

= (1| {p-^6+ (l2„, + 4iZ)P6'} ® |p-^6+ (l2„ + 4iZ)P6'} |NESS) 

= P-^(l|6 ® fe'lNESS) [(l2„ + 4iZ)P]^ = l2„ + 4iZ (93) 

where we used ([75]), and (1|6 ® S'|NESS) = l2n due to almost-CAR (^B>- | 

5. Examples 

We conclude this note by showing in two physical examples that the situations of non- 
diagonalizable Liouvilleans and subspaces of degenerate NESS are not as pathological 
or difficult to encounter as it might seem from the first sight. 

5.1. A single qubit with non-diagonalizable Liouvillean dynamics 

Take a single fermion (or qubit), described in terms of a pair of Hermitian anticomuting 
operators Wi = c + c^ and W2 = i(c — c^). Take Hamiltonian H = ihwiW2 and a single 
Lindblad operator L = vT(wi + e^^W2) Then, the key matrices read 

.. ^f r Tcos9 + h\ ,, / -Tsm9\ 

For h = Fcos^ the matrix X becomes nondiagonalizable. 

5.2. A pair of spins with degenerate NESS 

Consider a pair {n = 2) of qubits, or spins 1/2, described by Pauli operators ai 21 
coupled via Ising interaction H = }^Ja\a\ and attached to a Lindblad reservoir only 
through one of the qubits Li = Via^ ,L2 = r2crf . This model is fermionized in terms 
of Majorana operators Wi,2,3,4 P] using Jordan- Wigner transformation H = —^Jw2W3 
and the key matrices read 



X 









o\ 
-J 



\ 


J 




0/ 



M,^^ 



/ 1 





^ 


-1 

















\ 





0/ 



(95) 



where r± = r2 ± Fi. The rapidity spectrum is 

/3i = r+, (3,^, = ^±Jf^^ -.P, (3, = 0. (96) 

The Lyapunov equation (136!) for the steady state 2-point correlator (pO!) has a unique 
solution, in spite of zero eigenvalue (since it has multiplicity 1, see the proof of lemma 
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namely 

/ r+ J 0\ 

o"P "P r\ r\ r\ 

(97) 





/ 


r+ 


J 


°^ 


2r_ 


-r+ 
-J 














2r^ + J2 




\ 








0/ 



However, the non-uniqueness of |NESS) is manifested in non- vanishing and non- 
unique 1— point functions, namely since the non-unique NESS is parametrized as 
|NESS;a) = |NESS) + a|NESS; 4, 1), where |NESS;4,1) = 6^|NESS) = cl|NESS), 
we have {wj) = dj^^a. Using the (generalized) Wick theorem, we see that even 
order monomials are insensitive to non-uniqueness parameter a, whereas odd order 
monomials depend on a. The range of admissible values of a is obtained for 
example by exphcit construction of the spectrum of the NESS density matrix Pness = 
4X^7 ejoi} (l|a'i\^?2^T^3^i*4|NESS; a;)w7^u^2^u^?^4* ^^d require it to lie in the interval 

[0,1]: 

Richer structure of degenerate NESS is obtained by adding further spins 1/2 to such 
a one-sided open Ising chain, namely each new spin coupled to the rest of the chain via 
Ising interaction adds a new pair of completely imaginary structure matrix eigenvalues 
±ib (plus a pair of fully complex eigenvalues). 
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